Spin-isospin excitations and (3^ / EC half-lives of medium-mass deformed nuclei 

P. Sarriguren, E. Moya de Guerra, A. Escuderos 
Instituto de Estructura de la Materia, Consejo Superior de InvesUgaciones Cientificas, Serrano 123, E-28006 Madrid, Spain 

A selfconsistent approach based on a deformed HF+BCS+QRPA method with density-dependent 
Skyrme forces is used to describe f5'^-Aec&y properties in even-even deformed proton rich nuclei. 
Residual spin-isospin forces are included in the particle-hole and particle-particle channels. The 
quasiparticle basis contains neutron-neutron and proton-proton pairing correlations in the BCS 
approach, while neutron-proton pairing interaction is treated as a residual force in QRPA. We 
discuss the sensitivity of Gamow- Teller strength distributions and (3'^ /EC half-lives to deformation, 
pairing and the strength of the particle-particle interaction. The dependence on deformation is also 
compared to that of spin Ml strength distributions. 

O ' I- INTRODUCTION 

o 

. Exploring the properties of nuclei with unusual proton to neutron ratios is nowadays one of the major tasks and 
most active topics in nuclear structure physics, both theoretical and experimentally In this context, the beta- 
^ ■ decay properties of these exotic objects are of extreme importance in two different aspects. First, most nuclei of 

■ astrophysical interest [|| are those far from stability and their beta-decay rates have to be estimated theoretically. 
Second, beta-decay properties and nuclear structure are intimately related. It is clear that a precise and reliable 
description of the ground state of the parent nucleus and of the states populated in the daughter nucleus is necessary 
to obtain a good description of the /J— strength distribution, and vice versa, failures to describe such distributions 
would indicate that an improvement of the theoretical formalism is needed. 

^Jq ' Microscopic models to describe the /^-strength based on spherical single-particle wave functions and energies with 
ps) , pairing and residual interactions treated in Random Phase Approximation (RPA) or Quasiparticle Random Phase 
• Approximation (QRPA) were first studied in Ref. Extensions of these models to deal with deformed nuclei 
were done in Ref. [Q, where a Nilsson potential was used to generate single-particle orbitals. Subsequent extensions 

■ including Woods-Saxon type potentials [|| , residual interactions in the particle-particle channel |^ , consistent Hartree- 
Fock (HF) mean field with residual interactions treated in Tamm Dancoff approximation selfconsistent approaches 
in spherical neutron-rich nuclei ^ and based on an energy-density functional can be also found in the literature. 

1^ ' In a previous work we studied ground state and /3-decay properties of exotic nuclei on the basis of a deformed 

I , selfconsistent HF-I-BCS-I-QRPA calculation with density dependent effective interactions of Skyrme type. This is a 
^ ' well founded approach that has been very successful in the description of spherical and deformed nuclei within the 
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^ valley of stability 1 12 1 . In this method once the parameters of the effective Skyrme interaction are determined, basically 
^ I by fits to global properties in spherical nuclei over the nuclear chart, and the gap parameters of the pairing force 
are specified, there are no free parameters left. Both the residual interaction and the mean field are consistently 
• i-H obtained from the same 2-body force. This is therefore a reliable method, suitable for extrapolations into the unstable 
^\ regions approaching the drip lines. It is worth investigating whether these powerful tools designed to account for the 
properties of stable nuclei are still valid when approaching the drip lines. 

The residual interaction introduced in Refs. ||l^,|ll|, consistent with the mean field, acts in the particle-hole channel 
{ph). However, it has been pointed out (see for instance Ref. |]l3| and refs. therein) that for a complete description 
of the Z?"*" and /3/3 strengths, the inclusion of the particle-particle {pp) residual interaction is required. We gave in 
an example of the sensitivity of the results to the inclusion of such a force but a complete calculation within our 
formalism including this force is still missing. One of the purposes of this work is to extend our previous calculations 
by including a residual separable force in the pp channel and to study its effect in a systematic way to distinguish 
what is general and what is particular in the behavior of the nuclear response. To this end we include in the QRPA 
calculation a separable neutron-proton pairing force as part of the residual interaction. 

Following the same criteria as in our previous work |ll[| , we apply this formalism to the study of a series of proton 
rich isotopes approaching N = Z 'm the mass region A = 70. There are several reasons why this is a region of special 
interest to study /3-decay. One is that the Q-value of the decay {Qec) is quite large in proton rich nuclei ||l^. This 
means that a large fraction of the Gamow- Teller (GT) strength is contained within the Qec window. Therefore, one 
can investigate directly the structure of the GT strength distribution by /3-decay measurements without dealing with 
other indirect methods to extract the GT strength such as (p, n) or {n,p) charge exchange reactions. Another aspect 
to stress is that the mass region A ~ 70 is at the border or beyond the scope of the full shell model calculations. 
Predictions for the strength distributions, half-lives, and other decay properties obtained from the present formalism 
are of special relevance in this mass region since they may be at the moment the most reliable calculations. The present 
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results can then be used as a guidance to further experimental searches. It is also worth mentioning that this mass 
region is known to present particular structure effects, such as shape coexistence, which are the result of a delicate 
balance of competing nuclear shapes . It is then important to explore whether the dependence of the beta-decay 
properties on the deformation of the parent nucleus could be used as an additional piece of information to elucidate the 
shape of the nucleus. As a last point, we mention that, by approaching systematically the N — Z isotopes in various 
isotopic chains, we are placed in optimum conditions to observe whether agreement with experiment breaks down 
at some point as we approach N — Z . These isotopic chains are a laboratory to test the validity of our formalism 
and to look for failures related to the appearance of new physics such as neutron-proton pairing correlations not 
taken into account explicitly in the present mean field calculations, that may be important in N = Z nuclei. In this 
work we extend the study of Refs. ||l^,|ll| by including a residual neutron-proton pairing interaction in the J"^ = 1^ 
channel, by studying the influence of usual BCS pairing correlations in the N — Z isotopes, and by discussing the 
similarities between the GT strength distributions and the spin Ml strength distributions, which are the AT^ = 
isospin counterparts of the AT^ = ±1 GT transitions. 

In Section 2, we present a brief summary containing the basic points in our theoretical description. Section 3 contains 
the results obtained for the GT strength distributions with a discussion of the dependence on the residual interaction, 
pairing correlations, and deformation. We also compare our results to the experimental available information. Spin 
Ml strength distributions are studied in Section 4 discussing their analogies with the GT strength distributions. The 
conclusions are given in Section V. 



In this Section we summarize briefly the theory involved in the microscopic calculations presented in the next 
Sections. More details can be found in Refs. [ |lO|Jll| . Our method consists in a selfconsistent formalism based on 
a deformed Hartree-Fock mean field obtained with a Skyrme interaction including pairing correlations in the BCS 
approximation. We consider in this paper the force SG2 of Van Giai and Sagawa, that has been successfully tested 
against spin and isospin excitations in spherical and deformed nuclei p7| . Comparison to calculations obtained 
with other Skyrme forces have been made in Refs. flO| , p| , showing that the results do not differ in a significant way. 
The single particle energies, wave functions, and occupations probabilities are generated from this mean field. 

For the solution of the HF equations we follow the McMaster procedure that is based in the formalism developed in 
Ref. jl^ as described in Ref . . Time reversal and axial symmetry are assumed. The single-particle wave functions 
are expanded in terms of the eigenstates of an axially symmetric harmonic oscillator in cylindrical coordinates. We 
use eleven major shells. The method also includes pairing between like nucleons in the BCS approximation with fixed 
gap parameters for protons A,r, and neutrons A^, which are determined phenomenologically from the odd-even mass 
differences through a symmetric five term formula involving the experimental binding energies pO| ] . The values used 
in this work are the same as those given in Ref. . 

In a previous work pT[ we analyzed the energy surfaces as a function of deformation for all the isotopes under 
study here. For that purpose, we performed constrained HF calculations with a quadrupole constraint pl| ] and we 
minimized the HF energy under the constraint of keeping fixed the nuclear deformation. Calculations in this paper 
are performed for the equilibrium shapes of each nucleus, that is, for the solutions, in general deformed, for which we 
obtained minima in the energy surfaces. Most of these nuclei present oblate and prolate equilibrium shapes that 
are very close in energy. 

We add to the mean field a spin-isospin residual interaction, which is expected to be the most important residual 
interaction to describe GT transitions. This interaction contains two parts. The particle-hole part is responsible for 
the position and structure of the GT resonance HJllI and is derived selfconsistently from the same energy density 
functional (and Skyrme interaction) as the HF equation, in terms of the second derivatives of the energy density 
functional with respect to the one-body densities p^ ]. The residual interaction is finally written in a separable form 
by averaging the Landau-Migdal resulting force over the nuclear volume. 



II. THEORETICAL APPROACH 
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The coupling strength is given by (f^JT^ 
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and it is completely determined by the Skyrme parameters to,ti,t2,t3, a, the nuclear radius R, and the Fermi mo- 
mentum kp, obtained with the same Skyrme force for the given nucleus. 

The particle-particle part is a neutron-proton pairing force in the J'^ = 1"'' coupling channel. We introduce this 
interaction in the usual way P, p^ , p3| , that is, in terms of a separable force with a coupling constant n^j^, which is 
fitted to the phenomenology, 
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where 
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The two forces ph and pp in Eqs. (2.1) and (2.4) are defined with a positive and a negative sign, respectively, 
according to their repulsive and attractive character, so that the coupling strengths x k take positive values. 

The peak of the GT resonance is almost insensitive to the pp force and k^j, is usually adjusted to reproduce 
the half-lives However, one should be careful with the choice of this coupling constant. Since the pp force is 
introduced independently of the mean field, if k^j, is strong enough it may happen that the QRPA collapses, because 
the condition that the ground state be stable against the corresponding mode is not fulfilled. This happens because 
the pp force, being an attractive force, makes the GT strength to be pushed down to lower energies with increasing 
values of k^x- 

A careful search of the optimal strength can certainly be done for each particular case, but this is not our purpose 
here. Instead, we have chosen the same coupling constant for all nuclei considered here. This value has been obtained 
under the requirements of improving in general the agreement with the experimental half-lives while being still far 
from the values leading to the collapse. This will be discussed in more detail in the next Section. 

The proton-neutron quasiparticle random phase approximation (QRPA) phonon operator for GT excitations in 
even-even nuclei is written as 
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where a"*" {a) are quasiparticle creation (annihilation) operators, ljk are the excitation energies, and X'^^ , Y^jf the 
forward and backward amplitudes, respectively. The equations of motion are solved in the proton-neutron QRPA. 
From the QRPA equations the forward and backward amplitudes are obtained as 
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with Etyu — Ev + the two-quasiparticle excitation energy in terms of the quasiparticle energies Ei. M"^^ are given 
by 



(2.9) 
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with 
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where v's are occupation ampHtudes {u^ = 1 — v^) and spin matrix elements connecting neutron and proton 
states with spin operators 
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The solution of the QRPA equations can be found solving first a dispersion relation, which is now of fourth order 
in the excitation energies lo. Then, for each value of the energy the amplitudes M are determined by using the 
normalization condition of the phonon amplitudes 
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The technical procedure to solve these QRPA equations is well described in Ref. [P3| . 

In the intrinsic frame the GT transition amplitudes connecting the QRPA ground state |0) (Xujk |0) — 0) to one 
phonon states \lok) {j^tiK 1^) ~ \^k)), are given by 



{ujK\<JKt^\Q) = TM'^' 
The Ikeda sum rule is always fulfilled in our calculations 
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III. DECAY PROPERTIES 

The bulk properties of the isotope chains considered here were already studied in Ref. [ pT| . Binding energies, 
charge radii, quadrupole moments, and moments of inertia were discussed and compared successfully to available 
experimental data. Here we will concentrate on the decay properties. 

In this Section we discuss the results obtained for the GT strength distributions, half-lives, and summed strengths 
in the proton rich isotopes belonging to the ^ ~ 70 mass region. The results correspond to QRPA calculations with 
the Skyrme force SG2 and they have been performed for the nuclear shapes that minimize the HF energy. 

Figures showing the GT strength distributions are plotted versus the excitation energy of daughter nucleus. The 
distributions of the GT strength have been folded with F = 1 MeV width Gaussians to facilitate the comparison 
among the various calculations, so that the original discrete spectrum is transformed into a continuous profile. These 
distributions are given in units of g^/Air and one should keep in mind that a quenching of the gA factor, typically 
9A,eff = (0.7 — 0.8) g A, free IS cxpccted on the basis of the observed quenching in charge exchange reactions and spin 
Ml transitions in stable nuclei, where gs,eff is also known to be approximately 0.7 gsjree- Therefore, a reduction 
factor of about two is expected in these strength distributions in order to compare with experiment. This factor is of 
course taken into account when comparing to experimental half-lives. 



A. The role of the residual interactions 



Figure 1 illustrates the effect of the residual interactions on the uncorrelated 2-quasiparticle calculation. The 
calculations are done for the oblate and prolate shapes of ^"^Kr. The coupling strength of the ph residual interaction 
is obtained from Eq. (2.3), and its value for A = 7A and Skyrme force SG2 is Xgt' ~ ^■^'^ MeV. The coupling 
strength of the pp residual interaction is varied from K^lp = to Kqj, = 0.07 MeV. 

If we compare first the calculation with only the ph residual interaction (dotted line) to the uncorrelated one (thin 
solid fine), we can see that the repulsive ph force introduces two types of effects. One is a shift of the GT strength 
to higher excitation energies with the corresponding displacement of the position of the GT resonance. The other is 
a reduction of the total GT strength. Obviously these effects are more pronounced as we increase the value of the 
coupling strength Xgt- ^'^'^ introduce a pp residual interaction (dashed and solid lines), we can see that its 

effect, being an attractive force, is to shift the strength to lower excitation energies, reducing the total GT strength 
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as well. The shift and reduction of strength is more pronounced at high excitation energies and the position of the 
GT resonance is hardly affected by this interaction. The GT strength is also pushed a bit by the pp interaction to 
lower energies in the low energy excitation region. This effect, although small, is of great relevance in the calculation 
of the (3'^ /EC half- lives, which are only sensitive to the distribution of the strength contained in the energy region 
below the Qec window. By comparing the curves obtained from K^rp — 0.05 MeV (dashed line) and n^lp — 0.07 MeV 
(solid line), we can also see how these effects are more pronounced as we increase the value of the coupling strength. 
This is the reason why the usual procedure to fit these two coupling strengths Xgt ^^'^ 

follows: First one 

chooses Xgt to reproduce the position of the GT resonance, usually determined from charge exchange reactions {p, n) 
and {n,p), and then one chooses k^j, to reproduce the half-hfe. In our work, since the value of Xgt' determined 
selfconsistently and the known GT resonances are reasonably well described |Tl|] , we do not carry out this case by case 
fitting procedure. However, to determine the value of k^j,, we calculate first GT strength distributions and half- lives 
and compare the latter with experiment to extract a value that producing a reasonable agreement is still within the 
range of values compatible with the correct treatment of the QRPA. 

In Figure 2 we can see the result of the calculation of the half- lives in ^^Sr as a function of the coupling strength 
K^j,. The half-lives are obtained from the familiar expression, 

= ^ E / ^) 1(1- II II 0^) I' ' (3.1) 

where / (Z, m) are the Fermi integrals. Note that in our previous work [ pT[ , we used for these integrals the values 
tabulated in |24[| . In the present calculations we compute the Fermi integrals numerically for each Z,uj values. By 
this procedure we get more accurate results. 

We use D = 6200 s and include standard effective factors 
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i9A/gv)eff = 0-7n9A/9v)free =0.90. (3.2) 



The half-life decreases with increasing values of n^rp. This is clear because as k^j, increases the strength becomes 
more concentrated at low excitation energies below Qec ^md therefore the half- lives are smaller. This is true up to 
values around K^'rp — 0.1 MeV, where the collapse of the QRPA takes place. In the case of ^^Sr, we get an optimum 
value of K^'rp = 0.03 MeV and — 0.07 MeV for the oblate and prolate shapes, respectively. This value will 
of course depend, among other factors, on the nucleus, shape, and Skyrme interaction and a case by case fitting 
procedure could be carried out. Nevertheless, we have done calculations for other cases and found that, in general, 
values around k^Ij, = 0.07 MeV improve the agreement with experiment in most cases and since this is a valid value 
far from collapse, we have chosen k^^j, — 0.07 MeV as the value of the coupling constant of the pp residual interaction. 

In Ref. 1^ Homma and collaborators studied /3— decay properties using Nilsson-|-BCS-|-QRPA with ph and pp 
separable residual interactions. They considered a wide range of nuclei to extract phenomenologically the coupling 
strength of such residual interactions by fits to experimental half-lives. They obtained the following dependence on 
the mass number: Xgt = 5.2/A° '^ MeV, k^j, = 0.58/^° '^ MeV. For A = 70 this gives a ph strength Xgt = 0-27 
MeV and a pp strength Kq'j, = 0.03 MeV. The ph strength is somewhat smaller than the one obtained from Eq. (^.3|), 
derived selfconsistently from our mean field. The pp strength is also smaller than the value adopted here. This is 
consistent with the trend observed in Fig. 1 of Homma et al. ^ that shows that for decreasing values of x^t '-'^^ 
needs smaller values of Kq'j'. On the other hand, a discussion of the adequacy of our ph strength is also demonstrated 
in Ref. , where we compared our results on the position of the GT resonance with experimental data from {p, n) 
and (n,p) reactions in this mass region. The fact that our ph and pp strengths are somewhat different from those in 
Ref. [|j is not surprising since the mean fields are also different. 



B. The role of BCS correlations 



As already mentioned in the Introduction, our theoretical treatment does not explicitly include neutron-proton 
pairing correlations in the mean field. Thus, our quasiparticle basis only includes neutron-neutron and proton-proton 
pairing correlations in the BCS approach. In principle, one could extend the BGS treatment to include also neutron- 
proton pairing correlations in the mean field. This may be important particularly for N=Z nuclei. 

In Ref. PH ] we studied bulk properties (binding energies, charge radii, quadrupole moments, moments of inertia,...) 
of the nuclei considered here and found that the agreement between theory and experiment is as good for the N = Z 
as for the N = Z + 2, Z + A, Z + 6 isotopes. We concluded that the effect of neutron-proton pairing correlations on 
these bulk properties is roughly taken into account by the use of the phenomenological gap parameters A^r, A^. 
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The main effect of taking into account neutron-proton pairing in the HF+BCS calculation would be to increase the 
diffuseness of the Fermi surface [psf . This diffuseness is proportional to the gap parameters. It is therefore interesting 
to study the sensitivity of the GT strength to the gap parameters in the N = Z nuclei ^"^Ge, ^*Se, ''^Kr, and ''^Sr. To 
this end we compare in Fig. 3 the QRPA results obtained with the SG2 force in those nuclei for various values of the 
gap parameters differing by ±0.5 MeV from the values extracted from the phenomenology. To make the discussion 
easier we did not include the pp residual interaction in those figures. 

The main effect of the BCS correlations is to create new transitions that were forbidden in the absence of such 
correlations. Since the occupation probabilities are now different from or 1, the already existing peaks at A = 
will decrease when A > 0, while new strength will appear at other energies and will increase with increasing gap 
parameters. This new strength is in general located at high excitation energy, while the strength already present at 
A = is mainly concentrated at lower energy. As a consequence, the main effect of increasing the Fermi diffuseness 
is to smooth out the profile of the GT strength distribution, increasing the strength at high energies and decreasing 
the strength at low energies. 

It is clear from Fig. 3 that each particular case has its own peculiarities but the general trend described above plays 
the major role. It is also interesting to remark that in the high energy tail beyond the GT resonance, the effects of 
increasing A and increasing kF^j, are opposite. As we mentioned before, increasing nF^rp tends to deplete those tails. 

C. The role of deformation 

The role of deformation on the GT strength distributions can be summarized, in general, in two types of effects. 
First, deformation breaks down the degeneracy of the spherical shells and this implies that the GT strength distri- 
butions corresponding to deformed nuclear shapes will be much more fragmented than the corresponding spherical 
ones. Second, the energy levels of deformed orbitals coming from different spherical shells cross each other in a way 
that depends on the magnitude of the quadrupole deformation as well as on the oblate or prolate character. This 
level crossing may lead in some instances to similar profiles in the GT strength distributions of the various coexisting 
nuclear shapes but in other cases it may lead to sizeable differences between the GT strength distributions corre- 
sponding to different shapes of the same parent nucleus. This fact can be exploited to gain information on what is 
the nuclear shape of a nucleus by just looking at the structure of its /3-decay. 

Fig. 4 summarizes the main results of this work. They correspond to the GT strength distributions obtained from 
our HF-l-BCS-f QRPA with SG2 in the isotope chains of Ge, Se, Kr, and Sr. We can see in this figure the results 
obtained for the possible nuclear shapes of the isotopes ranging from N = Z to N = Z + 6. As usual the strengths 
are in units of g^/An and are plotted versus the excitation energy of the daughter nucleus. 

We can see that in general the N — Z isotopes of each chain contain the maximum strength as it corresponds to the 
most unstable nuclei. This strength becomes smaller and smaller as we increase the number of neutrons approaching 
the stable isotopes. The excitation energy of the GT resonance also decreases with the number of neutrons. 

If we look in Fig. 3 the case of Ge and Se isotopes, we can see that there are no signatures of the actual character, 
oblate or prolate, of the deformation of the parent nucleus: Both equilibrium shapes produce similar GT profiles 
independently of whether it is oblate or prolate. Then, it can be concluded that in these isotopes one cannot use GT 
distributions to distinguish between the two coexistent shapes predicted by the theory. On the other hand, the figures 
corresponding to Kr and Sr isotopes, show differences in the GT profiles of the various shapes for each isotope that 
can be easily identified even within the Qec window. A case by case analysis allows one to conclude that the most 
favorable candidates to look for deformation effects based on /3-decay measurements are '^^Ki and '''6.78,80^^.^ these 
isotopes the GT strengths within the Qec window are different enough to distinguish between different equilibrium 
shapes. On the contrary, there are cases like ^^Kr where though the large Qec value {Qec = 5 MeV) make it 
worth to investigate experimentally, the measured strength may not be conclusive on the shape because oblate and 
prolate shapes produce similar profiles. Other isotopes like 76.78]^j. j^^ve small Qec values that will not allow a clear 
conclusion either, while ^^Kr seems to be the best Kr candidate. It has a Qec window (3.1 MeV) big enough to 
distinguish between a continuously increasing profile, as the prolate shape predicts, or a completely developed bump 
structure peaked at around 1.5 MeV and vanishing at about the Qec value, as the oblate shape predicts. A similar 
situation happens in the case of ^^Sr {Qec — 6.1 MeV). The oblate shape produces a peak centered at around 1 
MeV and a second one centered at 6 MeV, while the prolate shape produces a GT distribution that increases almost 
continuously up to 6 MeV. The case of ^^Sr {Qec — 3.8 MeV) is again a clear case where one can distinguish between 
the bump structure generated by the spherical shape or the continuously increasing pattern of the prolate case. The 
same is true for the *°Sr isotope, although in this case Qec is not so big {Qec = 1-9 MeV). 
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D. Comparison to experiment 



In this subsection we compare to experiment our calculated Qec a-nd total /EC half- lives Ti^2y well as a few 
cases where the GT strength distribution have been partially measured by /3-decay. We also discuss the GT strengths 
contained within the Qec windows. The difference with respect to the values quoted in our previous work is that 
now the pp interaction is included. There is also a minor difference in the calculation of Qec, 

Qec = ™7r — + — {\u + E^)^j^_^^ + (A„ — £'„)(jv,z-2) ■ (3-3) 

The Q EC values in Tables 1-4 of Ref. were calculated approximating the lowest quasiparticle energies E^, and E^^ 
by the gap parameters A^, and Att, respectively, while Qec values in Table 1 here are calculated using their actual 

values E = ^J{e-\f + 

We can see in Table 1 these results obtained from QRPA calculations with the Skyrme force SG2 and for the different 
shapes, oblate, prolate or spherical, where the minima occur for each isotope. We did not include the stable nuclei, 
all having Qec less than zero and infinite half- lives. It should also be mentioned that the total GT strength below 
Qec has been reduced by a quenching factor 0.6 to be consistent with the same quenching used in the calculation of 



the half- lives in Eq. (3.1) 



We get a very good agreement with the measured Q ec values in practically all cases and the values obtained with 
the various shapes are quite similar. The half-lives are also well reproduced in most cases. Only in the most stable 
isotopes, where the half-lives are very large, we find some noticeable discrepancies but this is not so relevant because 
in these cases the Qec values are very small and therefore, the half- lives are only sensitive to a tiny region of the GT 
tail at low energies. In other cases the agreement is very reasonable. The sums of the GT strength up to the Qec 
value do not differ much from one shape to another although this does not mean that the structure of the profiles are 
equivalent. As we have seen in the last subsection, there are cases, ^^Kr and 76,78,80gj.^ whose profiles can be easily 
distinguished although the final summed strength is very similar. 



Table 2 contains experimental information on GT summed strengths in the N — Z nuclei Kr pq] and Sr |27| up 



to different values of the excitation energy, always below Qec- They are compared to our theoretical calculations for 



the two equilibrium shapes using the same quenching factor as in Eq. (3.2). From this comparison an oblate shape 
for ^^Kr and a prolate shape for ^^Sr seems to be favored, but this is not yet sufficient for a conclusive answer. 

In any case it would be worth to extend the measurements in ^^Sr to higher excitation energies because it is a very 
good example where deformation effects are visible (see Fig. 4 for ^^Sr). On the other hand, the extension of these 
measurements to the case of ^^Kr may not be so relevant since one could not distinguish one shape from another (see 
Fig. 4 for 72Kr). 



IV. MAGNETIC PROPERTIES 



The spin Ml transitions are the AT^ = isospin counterparts of the AT^ = ±1 GT transitions. The study carried 
out in |1^] for several isotope chains showed that the Ml strength depends on deformation, not only for orbital but 
also for spin excitations. Thus, on general grounds, one may argue that the deformation dependence of the GT 
strength discussed earlier should be comparable to that of the spin Ml strength considered in |1^. This is expected 
to be particularly the case when N = Z. 

In a similar way to Fig. 4 for the GT distributions, we show in Fig. 5 the profiles of the spin Ml strength 
distributions for the same isotope chains. The results correspond to the selfconsistent HF-I-BCS-I-QRPA calculations 
with the SG2 interaction. The details of these calculations are described in and closely follow those of the GT 
strength except for the ATz = character of the Ml operator. 

Similarly to the GT strength distributions in the Ge and Se isotope chains in Fig. 4, we can see that the Bcr{Ml) 
strength distributions for these two isotopes in Fig. 5 have similar structures for the oblate and prolate shapes. They 
have a big resonance located at the same energy around 10 MeV and a small bump at about 5 MeV. However, the 
strength contained in the prolate peak is always larger than the oblate one, contrary to what happened with the GT 
strengths that were comparable. The spherical shape in ^^Se produces much less strength than the deformed shapes. 
We can also see in Fig. 5 for Kr and Sr isotopes that the profiles of the Ml strengths corresponding to the prolate 
and oblate distributions can be clearly distinguished, similarly to the case of the GT distributions in Fig. 4. The 
strength corresponding to the prolate shape is again the largest. Therefore, clear similarities between the GT and 
Ml strength distributions are observed, but some differences can be seen as well. In particular, for spin Ml strength 
distributions the position and strength of the resonance are practically the same in all nuclei in a given isotope chain. 
This is different to what happened with the GT strength distributions, where the /3+-GT strength decreases very fast 
with increasing N — Z due to Pauli blocking, while Ml strengths are not affected by this. 
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V. CONCLUSIONS 



We studied /S'^-decay in various isotopic chains of medium-mass proton-rich nuclei within the framework of the 
selfconsistent deformed HF-I-BCS-I-QRPA with Skyrme interactions. Our spin-isospin residual interaction contains 
a particle-hole part, which is derived selfconsistently from the Skyrme force, and a particle-particle part, which is 
a separable force representing a neutron-proton pairing force. Compared to previous calculations, where the latter 
residual interaction was not considered, we obtain in general better agreement with experiment for half-lives. It is 
worth mentioning that for the N=Z isotopes this agreement is comparable to that achieved for the other isotopes 
with an excess of neutrons. This indicates that using the phenomenological gap values for neutrons and protons, as 
well as the neutron-proton pairing correlations as a residual force in QRPA, is sufficient to account at least for this 
experimental information. 

From our study of the dependence on the shape of the GT strength distributions we conclude that there are some 
interesting cases that are worth to explore experimentally such as ^^Kr and 76,78,80gj.^ these examples the measured 
P'^ strengths below Qec could be used to identify their equilibrium shapes. In particular, we have observed in ''^Sr 
various compatible indications pointing out towards the same conclusion: a strong prolate component in the ground 
state that is suggested from the comparison to experiment of half-lives and GT strength measured at low excitation 
energies. Nevertheless, more experimental information is still needed to reach a conclusive answer. 

We also studied the energy distributions of the spin Ml strength and the similarities and differences with their GT 
counterparts. We conclude that the main features of GT and spin M 1 strength distributions are similar. This suggests 
that we may learn about properties that are observable in highly unstable nuclei (like (3 decay) from properties that 
are observable in stable nuclei (like Mi's), and vice versa. Information on spin Ml excitations can be extracted from 
both leptonic and hadronic probes. In stable nuclei, the combined analysis of electron, photon, and proton scattering 
data has provided reliable information on orbital and spin Ml strength distributions ||2^. In principle, the same type 
of experiments could be performed on the proton rich nuclei considered here. Inelastic {p,p') scattering experiments 
might be particularly suitable to explore the spin part of the Ml strength and to test the predictions of the present 
paper. 
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Table 1. Results for Qec values, half-lives {T1/2), and GT strength summed up to Qec energies (YIiq^^)- The results 
correspond to QRPA calculations performed with the Skyrme force SG2 in four isotopic chains and are calculated for 
the various equilibrium shapes in each nucleus. Experimental values for Qec and T112 are from pO[ |. 







Opt 




T 

1/2 




V B(GT) 






exp 


th 


exp 


th 




ODl 


/I /1 1 




00. ( b 


o4.0 s 


n 7 




pro 




A 1 




ID ( .U S 


u.o 




obi 


2.10 


2.3 


2.3 h 


1.6 h 


0.3 




pro 




9 9 




1 Vi 

0. -L 11 


n 9 

u.z 


68 Qg 


obi 


0.11 


0.2 


271 d 


198 d 


0.1 




pro 




u.o 




lUU 11 


n n 

u.u 




r»hl 
UUi 








77 9 s 


1 1 

L . L 




pro 




4.0 




DD.4 S 


n Q 
u.y 


obi 


2.40 


2.5 


41.1 m 


38.8 m 


0.5 




pro 




9 7 

z. f 




OrJ.iJ 111 


u.o 








u.o 


S 4 H 


H 

0.0 u 


n 1 

U. 1 




pro 




1 

1.0 




n H 

u.o u 


n 9 

u.z 


XVi 


(JUi 






1 7 9 e 


91 4 s 

Z 1 .4; b 


1 9 
1 .z 




pro 




9 




10. D S 


1 Q 




obi 


3.14 


3.4 


11.5 m 


8.7 m 


0.7 




pro 




3.5 




12.4 m 


0.6 


^«Kr 


sph 


1.31 


1.7 


14.8 h 


4.1 h 


0.2 




pro 




1.2 




38.0 h 


0.1 




obi 


6.10 


5.9 


8.9 s 


3.2 b 


2.1 




pro 




5.8 




10.9 b 


2.3 


^«Sr 


sph 


3.76 


4.3 


2.7 m 


1.3 m 


0.4 




pro 




3.1 




19.9 m 


0.6 


«OSr 


sph 


1.87 


1.8 


1.8 h 


56.0 h 


0.1 




pro 




1.6 




6.4 h 


0.2 



15 



Table 2. Comparison of the GT strengths contained below some given excitation energy between experimental 
measurements ( for ^^Kr, for ^^Sr) and theoretical calculations. 
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